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Abstract 

We re-examine approximations in the analytical calculation of the primor- 
dial spectrum of cosmological perturbation produced during inflation. Taking 
two inflation models (chaotic inflation and natural inflation) as examples, we 
numerically verify the accuracy of these approximations. 
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The inflationary prediction on the primordial power spectrum of the cosmological per- 
turbation is testable by the observations of the large-scale structures therefore it is 
important to make the prediction more precise. In this brief report, we re-examine the ana- 
lytical calculation of the primordial spectrum of cosmological perturbation produced during 
inflation and point out a subtlety in the approximations. 

To begin with, we review briefly the derivation of the power spectrum using the standard 
slow rolling approximation (SRA) technique [|2]||. 

The most general form of the metric for the background and the scalar metric perturba- 
tions can be written as [|J 

ds 2 = a 2 (r]){(l + A)drf - 2diBdx l dr] - [(1 - 2*)^- + 2d i d j E}dx i dx j } , (1) 

where r\ = / dt/a(t) is conformal time. 

The primordial spectrum of scalar perturbation is defined as 

p s (k) = ^- 2 \n k \ 2 , (2) 

where TZ k is the coefficient of the Fourier transformation of the intrinsic curvature pertur- 
bation of comoving hypersurfaces 1Z which during inflation has the form |3|,[|] 

K = -¥ - S-6<f> , (3) 

where 5<p represents the fluctuation of the inflaton field and H = a/ a is the Hubble 
constant. 

Note that Ps(k) in Eq.@ by definition is constant, however, as we can see later the 
power spectrum during the period of inflation evolves in time and approaches a constant 
only when k/aH — > 0. For the convenience of our discussion we introduce Ps{k,a), the 
power spectrum of a given k mode as a function of time, to distinguish from its asymptotic 
constant value P s (k) = P s (k,a)\ k / aH ^ . 

To calculate Ps{k, a), one usually defines u k = —zlZ k where z = a<p/H. And u k satisfies 
the following equation of motion || 
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4 + \ k --JUk = 0, (4) 

where a prime denotes the derivative with respect to conformal time rj. 

According to the flat spacetime field theory (well inside the horizon), one has 

u k -> -j^e~ ikr i as k/aH -> oo. (5) 

On the other hand, the solution of Eq.@ in the limit of k/aH — > is 

Uk oc z . (6) 

Eq.(|D indicates that Ps{k, a) becomes approximatly constant for mode k to be well outside 
the Hubble radius, which can also be seen directly from the equation of motion of TZf.: 

72* = !(l + «;)-W^ **, (7) 

here, w = p/p is the ratio of the pressure to the density. 

To describe the period of slow-rolling of inflation, one defines two parameters e and S 

-H 

< 1 , 



H 2 



z" 



(8) 

in term of which z"/z can be expressed as [BJ 

2a 2 ^ 2 (l + 16 + e + ±5 2 + \e5 + + & S) , (9) 
and the conformal time 

f dt f da —1 f eda 
r] = J^ = J^H = ^H + J^H- (10) 

Assuming e and 5 to be constant during the period of slow rolling, Eq. (|) can be reduced 
to I: 

Ufc+ffc 2 -^^)u* = , (11) 



where v = (1 + 5 + e)/(l - e) + 1/2 and 

We emphasize that to obtain Eqs.flllD and (|12"D, one has neglected the terms proportional 
to e or 5 in Eq.@ and taken e in Eq.fllOD constant. 
The solution of Eq.flTTD for iz fc is 

u k = l^ v+ ^(-ri^Hjp(-kri) . (13) 



In the limit of k/aH — > 0, becomes 



Substituting rj of Eq.(12) into Eq.(14), one has the power spectrum in the slow rolling 
approximation, 

^) (tal -^%- e r^W as ±^0 (15) 

To obtain the final value of power spectrum, what has been done in the literature is to 
evaluate the time dependent Pg (k,a) in the equation above at k — aH 0|§: 

P s SRA) (k)^Q(k,a)\ aH=k (16) 
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47T 2 > 

where Q(k,a) is defined as the right hand side (R.H.S.) of Eq. flT5|) . 

Note that P { s SRA) (k, a) in Eq.flTo]) is valid only in the limit of k/aH — > 0, however, 
to obtain Eq.(|T7D one takes k/aH = 1. One wonders if the value in Eq.flTTD is a good 
approximation to the asymptotic value Ps{k) in the limit of k/aH — > 0? Furthermore, 
<5(&, a) is mathmetically divergent for z/ > | and zero for z/ < | in the limit of k/aH — > 0, 
so Q(/c, a) can only be a description of the true evolution of Ps{k, a) for some specific range 
of \n(aH/k). Therefore it is questionable to have the asymptotic value of the real spectrum 
obtained from Eq . (O) . 



In this brief report, we numerically solve Eq.(|j) and study the evolution of Q(k,a) and 
Ps{k, a) after horizon crossing. The initial conditions of Uk and u' k are given by Eq. ([I3|) when 
the mode k crosses out the Hubble radius. In the numerical calculation of the amplitude of 
the primordial spectrum we take two models as examples: one is chaotic inflation model with 
a potential V = m 2 2 /2, and the other is natural inflation model |7]] V = A 4 [l + cos(0//)]. 
We consider only the amplitude of the primodial spectrum for the COBE scale (kcoBE — 
7.5aoH ). The fluctuation mode corresponding to this scale crosses out the Hubble radius 
at the time e-folds iV ~ 57 before the end of inflation. This condition determines the value 
of inflaton field 4>\ a H=k CO BE- O ur numerical results on the evolution of Ps(kcoBE,a) and 
Q(k C oBE,a) are presented in Fig.l and Fig. 2. In each of these figures, a dashed horizontal 
line marks the exact asymptotic value of the power spectrum. 

Fig. 1(a) and Fig. 2(a) show the evolution of Ps{kcoBE, a ) after the COBE scale crosses 
outside the horizon during inflation. From these two figures one can see that when 
\n(aH/k C oBE) ~ 3, Ps(k C oBE, a ) starts to be approximatly a constant. In Fig.3 and 
Fig. 4 we show the relative error between QikcoBE,®) and the exact numerical result of 
Ps(kcoBE,a) when ln(aiT / kcoBE) > 3. For 3 < \n(aH/kcoBE) < 6 it is very small, so 
Q(k C oBE, a) in this range of \n(aH/kcoBE), which we denote as Q(k C oBE, a )l7t~o> wm ^ e a 
good approximation to the real spectrum Ps{kcoBE, a) in the same range, and then can be 
taken reasonablly as the asymptotic value of the spectrum in the limit of k/aH — > 0. 

Fig. 1(b) and Fig. 2(b) show the evolution of Q(kcoBE, a ) after crossing outside the hori- 
zon. From these figures, we can see that the uncertainty in calculating Q(k C oBE, a )\n^o ^ 
using the approximate formula Eq.(|17D is of order 1%. We conclude that the usual analytical 
treatment in obtaining the primordial power spectrum is feasible with a good accuracy. 

Before concluding this report, we should point out that Nakamura and Stewart || have 
tried to tackle this issue by expanding analytically the R.H.S. of Eq.(|15D at k = aH to the 
first-order term of ln(aH/k). However, as shown in Fig.l and Fig. 2, when Ps{kcoBE,a) 
starts to be approximately a constant, \n(aH/kcoBE) is larger than about 3. This raises 
another question about the validity of the analytical expansion, which can not be answered 
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by an analytical argument. So a numerical study presented in this paper is necessary. 
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FIG. 1. Evolution of Ps(kcoBE, «) and Q(kcoBE, a) as a function of In (aH/kcoBE) for model 
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FIG. 2. Evolution of Ps(kcOBE, «) and Q(kcoBE, a) as a function of In (aH/kcoBE) for model 
V = A 4 [l + cos(</>//)] with / = Imp,. 
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FIG. 3. The relative error between Q(kcoBE,a) and Ps(kcoBE,a) for model V(<j>) = m 2 (f> 2 /2 
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FIG. 4. The relative error between Q{kcoBE,a) and Ps(kcoBE,a) for model 

V = A 4 [l + cos(0//)] with / = lm p i . 
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